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Abstract. We derive discrete and oscillatory Chern-Simons matrix models. 
The method is based on fundamental properties of the associated orthogonal 
polynomials. As an application, we show that the discrete model allows to 
prove and extend the recently found relationship between Chern-Simons the- 
ory and g-deformed 2dYM. In addition, the equivalence of the Chern-Simons 
matrix models gives a complementary view on the equivalence of effective su- 
perpotentials in J\f = 1 gauge theories. 



1. Introduction 

Since the seventies, when random matrix theory [1] was successfully employed 
in gauge theory [2], matrix models have played a rather remarkable role in gauge 
theory (see [3, 4] for reviews). In this paper, we shall focus on the matrix models 
that appear in Chern-Simons theory [5, 6] and also discuss the recently found 
connection with 2dYM [7, 8] as well asAf—l supersymmetric gauge theories. More 
precisely, with recent work that discusses the equivalence of effective superpotentials 
in these theories [9, 10, 11]. 

Let us briefly recall very basic facts of Chern-Simons gauge theory: In [12], 
Witten considered a topological gauge theory for a connection on an arbitrary 
three-manifold M, based on the Chern-Simons action: 

(1) Scs(A) = -^ [ Tr(A A gL4 + \a A A A A), 

4tt J m 3 

with k an integer number. One of the most important aspects of Chern-Simons 
theory is that it provides a physical approach to three dimensional topology. In 
particular, it gives three-manifold invariants and knot invariants. For example, the 
partition function, 



(2) Z k {M) =1 VAe* Scs{A \ 

delivers a topological invariant of M, the so-called Reshetikhin-Turaev- Witten in- 
variant. Recent reviews are [4, 13]. 

As reviewed in detail in [13], a great deal of interest has been recently focused 
on the fact that Chern-Simons theory provides large N duals of topological strings. 
This connection between Chern-Simons theory and topological strings was already 
pointed out by Witten [14] (see also [15]), and then considerably extended in [16]. 

Interestingly enough, Marino has found a tight connection between Chern-Simons 
theory and random matrix models [5]. More precisely, building upon [17], in [5] it is 
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shown that the partition function of Chcrn-Simons theory on S 3 with gauge group 
U(N) is given by the partition function of the following random matrix model: 

(3) Z = 2^pl/ne-«^. n (2 Si nh^) 2 f, 

i=l i<j V 7 

From the point of view of topological strings, this describes open topological A 
strings on T* S 3 with N branes wrapping S 3 [5]. Chern-Simons matrix models 
have been further considered in [6] and [18]-[26]. Most of these works focus on the 
relevance to topological strings. In [6], the emphasis is on exact solutions and on the 
special features of the matrix models. In this paper we shall be further developing 
these results, while making some connections, at the end, with matrix models in 
M = 1 suspersymmetric gauge theories. 

It should be stressed that the result (3) in [5] is just a particular case of more 
general models and our results will apply to them as well. Explicit expressions for 
them will be given in the next section. 

As shown in [6], the Stieltjes-Wigert polynomials, a member of the q-deformed 
orthogonal polynomials family [27], allows to compute, in exact fashion, quantities 
associated to the matrix model. In the computation, the ^-parameter of the poly- 
nomials turns out to be naturally identified with the q-parameter of the quantum 
group invariants associated to the Chern-Simons theory. In [6], the focus was on 
U(N) Chern-Simons theory on S 3 , whose matrix model is also the Stieltjes-Wigert 
ensemble: 

(4) Z = J[dM]c-^t T ^ osM ^ . 

A special property of any of these type of matrix models is that there are infinitely 
many models that lead to the same partition function [6]. This is a courtesy of 
their very weakly confining potentials. We shall extend this result, essentially by 
exploiting results from the moment problem [28, 29], as in [6], in a more explicit 
way. In some cases, it just suffices to consider well-known facts of orthogonal 
polynomials with undetermined weight [29]. In others, the ones corresponding 
to the oscillatory/periodic models, one should further develop -in a simple and 
natural way-, the results of the moment problem to fully understand the extent of 
the equivalence between the models. 

To summarize, as a first result we have the existence of Chern-Simons matrix 
models with a discrete weight, which essentially turns out to be a discretization 
of the weight of the usual models [5, 6]. That is to say, the continuous models 
in [5, 6] can be castled into an homogenous and exponential lattice, respectively. 
As we shall see, this discrete model is precisely the q-deformed 2dYM partition 
function. The advantage of our method is that it can be applied to any Chern- 
Simons matrix model. Then, we shall deal with the continuous case, and show 
the manifold possibilities, including losing differentiability of the weight function. 
From a mathematical point of view, it is remarkable that the potential may be 
so wild even to reach fractal behavior and still leave the partition function (which 
in addition is a Chern-Simons partition function) unchanged. Interpretations in 
terms of gauge theories and quantum groups are discussed in some instances, and 
we show that the equivalence of potentials of the (generalized) Chern-Simons matrix 
models introduced, leads to a complementary view of the equivalence of effective 
superpotentials in N = 1 supersymmetric gauge theories [9, 10, 11]. In the last 
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section, we comment on the possible physical relevance and some further avenues 
of research are suggested. 

2. Generalized Chern-Simons matrix models 

The idea, already introduced in [6], is that Chern-Simons quantities can be ob- 
tained from a very precise coarse-grained knowledge of infinitely many different 
distributions, not just the log-normal. In matrix model language: there arc infin- 
itely many, equivalent, matrix models. Remarkable particular cases are: 

1) Discrete matrix model. The discretization of the models in [5, 6] is deeply 
related to the quantum group symmetry of Chern-Simons theory, and leads directly 
to 2dYM [7, 8]. 

2) The continuous and periodic models. The main feature is the presence of 
log-periodic [30, 31] 1 and periodic behavior. 

These examples underline to what extent the matrix models can be deformed 
while preserving the Chern-Simons information. 

2.1. Preliminaries. For future use, we quote here the moments of the weight 
function associated with the Stieltjes-Wigert polynomials: 

(5) 4= rdxx n e~ k2 ^ 2 =e (" +1 ) 2 / 4fe2 . 

Jo 

In what follows we will set q = e _1 / 2fe2 . Notice that the above formula is valid 
as long as Re/c 2 > 0. For the applications in Chern-Simons theory, we will be 
interested in the two cases Re k 2 = and Im k 2 = 0. 

We now consider the set of weights that satisfy the functional equation: 

(6) w{qx) = ^/qxw(x), 
it is easy to see that w also satisfies 

(7) w(xq n ) =q n2 / 2 xw{x). 
Let us normalize w such that 

/>oo 

(8) <T 1/4 / dxw(x) = l. 

Jo 

Applying the proof of proposition 2.1 of [35] to this weight, it now follows that it 
has the same moments as the log-normal weight: 

/>OC 

(9) / dxx n w(x) =<r («+i) 2 /2 . 

Jo 

The log-normal distribution (5) is a distinguished case of an undetermined moment 
problem and a basic result from the theory of moments is that the set of solutions of 
an indeterminate moment problem contains discrete measures as well as absolutely 
continuous measures. Indeed, it was Stieltjes in his seminal work [28], the first to 
show that all the functions in the following family: 

(10) f#(x) = e" los2;r (l + ^sin(27rloga;)), with i? e [-1, 1] , 

^This behavior is typical of models with a discrete-scale invariance [31] . It appears in rigor- 
ous formulations of fractal geometry [30] and in several condensed-matter physics applications 
[31]. Connections with quantum group symmetries are discussed in [32]. Recent physical appli- 
cations involving this behavior include limit-cycles in RG flows (see [33] for example) and critical 
phenomena in gravitational collapse [34]. 
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have the same moments, the log-normal moments. Thus, the parameter i? does 
not play any role regarding integer moments and the infinitely many functions in 
the family g$ (x) all have the same integer moments (and consequently, the same 
orthogonal polynomials). This case corresponds to k = 1. More generally, for an 
arbitrary k, the function is then: 



(11) 



ft(x) = e - fe2log2 * (l + i?sin^27r^!^ , with t? e [-1,1] 



and with q = e 1 / 2 ' c2 as above. 

We consider one more class of weight functions, namely the discrete ones [35, 36]: 

oo 

(12) W (x) = — — £ c"q^6(x of), 

V y ^ ' n= — oo 

where 

(13) M(c) = {-cq^/q, -y/q/c, q; q)^ . 

From now on we will take q to be real and c an arbitrary real positive constant, 
c > 0. It is not hard to see that w again has the same moments (9): one just uses 
the Jacobi triple product identity 

oc 

(14) V (-l) n q®x n = (x,q/x,q;q) 0O , 



n— — oo 

and translates the sum. 

2.2. Generic description. In fact, it is possible to have a more general charac- 
terization of the solutions of the log-normal moment problem. Following [36, 35], 
consider two generic weight functions, cj\ (x) and uj 2 (x) satisfying the functional 
equation (6). Since the functional equation determines the log-normal moments, 
these generic functions possess the log-normal moments. Consider the quotient: 
g (x) = Ldi (x) I L02 {x) , which leads to a function that satisfies a functional equa- 
tion: 

(15) g(x)=g(qx), 

and thus, two generic functions with log-normal moments differ by a function that 
satisfies (15). This type of function can be naturally named, as in [36, 35], q- 
periodic, and it is manifest that any logarithmic oscillatory term (with the proper 
frequency, as in (11) for example) satisfies this property. This type of function is 
very natural when q deformations are present. 

Thus, any integral constructed from the moments with any of the above weights 
will give the same result. We can now apply this to matrix models where we have 
N such integrals. Consider the following matrix model: 

/•CO />oo 

(16) Z= d Xl ... dx N e-^Z*=i lo z Xi l[(x k -x l ) 2 P(x u ...,x N ), 

J ° Jo k<i 

where P has a finite Laurent expansion in the x's. We write it as follows: 

(17) P{X U ...,X N )= ]T Pni ... nN xT ---X^, 

rii...n N 

for some coefficients p ni ... njv , where the n^'s can be positive or negative integers. 
Notice that P is a symmetric Laurent series in the a^'s, since an antisymmetric 
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piece in any two Xi's would not contribute to the integral. We can expand P in a 
basis of symmetric polynomials or, alternatively, a basis of Schur functions. 

It is now obvious that (16) inherits the invariance under deformation of the 
weights as discussed above. The Vandermonde interaction itself is a polynomial of 
degree N(N — 1) in x, so the integrand factorizes into a product of integrals over 
the Xj's, and the weights can be deformed for each integral. Thus, we get that 

/>oo />OC 

Z = / d Xl ... dx N e~?k^^ los x *l[(x k -x l ) 2 P(x 1 ,...,x N ) 
Jn Jo k<i 

(18) = / dxi... dx N wix^ . . .w(x N ) Y\(xk - xi) 2 P(x 1 , . . . ,x N ), 

Ja Ja k<l 

where w(x) is any of the above weights, whether continuous or discrete. 

Therefore, the following matrix model is also invariant under the same deforma- 
tions: 

(19) Z = y > [dM]e^537 Tr(logM)2 .f(M), 

where f(M) has a finite Laurent expansion in the moments TrM™, and M is a 
Hcrmitian matrix. 

We now apply the above to Chern-Simons theory. Consider the partition function 
of Chern-Simons theory on a Seifert space M ~ AT(^-, . . . , 2a.). This is obtained by 
doing surgery on a link in S 3 with n + 1 components, out of which n arc parallel, 
unlinked unknots, and one has link number 1 with each of the n unknots. The 
surgery data are Pj/qj for the unlinked unknots, j = 1, . . . , n, and for the last 
component. The partition function is 2 [5]: 

(20) Z CS( M) = (~1) |A+I [^P)f A ^ e ^ sienWP) -^ 
{ ' CS{ 1 \W\(2wi) r \ VolA ry l |P|<V 2 

V" J Mp -0 2 /2gs-lt-p n^iria>o 2si nh|^ 



d(3e 



teA t /HA t J nQ>o ( 2s inh^) ? " 



This expression gives the contribution of the reducible flat connections to the par- 
tition functions. Recall that for both S 3 and lens spaces this amounts to the exact 
partition function. The case n = corresponds to the three-sphere S 3 that leads 
to (3). Thus, for the case of U(N), and focusing on a particular sector of flat 
connections, we get the following matrix model: 

tV-cc ^(2 sinhM_m)" 

We used a different letter to remind ourselves that to obtain the full partition 
function one needs to include the constant prefactor and sum over t. Notice that 
because of the integrals, the integrand in the above expression is automatically 
symmetric in the tj's. 

All we have to do in order to show invariance under deformations of the weight 
is to show that this is of the form (16). But this should be clear. As in [6], we 
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perform a coordinate transformation 

(22) yi = a\ogx l + b, 

for some constants a and b. We get: 

2 t°° / i OC 2 

Zcs(M) = a N e-^t- bl ^ u dx x . . . dx N e~^^=^ Xi { Xl . . .x N )~ d 

Jo Jo 

( 23 ) >< n i n f (*2 - *?) a -" n(»2 /w - . 

i=l X i fc<2 V i=l / 

where 

7 -, TO , 

(24) d=l + - + ^a(N-l)(£ -+2-n). 

Obviously, a sufficient condition for this to be of the form (16) is to choose a = 
api . . .p n , b — cg si where a, b e Z. Thus, all these models can be deformed in 
the way shown above. Basically, any matrix model with an integrand where the 
periodicity is broken by a Gaussian weight, will have this deformation property. 
This includes all of the known Chern-Simons matrix models, the most general one 
being that for torus link invariants [37]. We will discuss this elsewhere [38]. 

2.3. Discrete model and 2d qYM. In the early 70's, Chihara employed the 
functional equation (6) to show the existence of an infinite number of discrete 
measures with the same log-normal moments [36]. An example is (12), discussed 
above. This case is remarkable, since it is completely explicit. Furthermore, it 
not only leads to a discrete Chern-Simons but it also turns out that this model is 
immediately relevant to the correspondence between Chern-Simons theory and 2D 
Yang-Mills [7] as we shall see at the end of this section. 

Note that the moments of (12) do not depend on the parameter c. Thus, as in 
the case of the continuous functions, this example already contains infinitely many 
functions with the same moments. The discrete case leads, for U(N) and 5 3 , to 
the following Chern-Simons matrix model : 

/N oo 
Hdxi cY +n2/ ^(s-«/ n )n^-^) 2 
i—l n— — oo i<j 

oo oo TV 

Cn £ ••• J2 f[c ni q n < +ni/2 l[(q n *-q n «) 2 . 

ni— — oo njy— — oo 2=1 i<k 

At this point, it is worth to recall the works on combinatorial quantization 
[39, 40], whose aim is a mathematically careful quantization of pure Chern-Simons 
theory. These works simulate Chern-Simons theory on a lattice, in such a way that 
partition functions and correlators of the lattice model coincide with those of the 
continuous model. Interestingly enough, in [39] it is shown that one can also choose 
the point of view of constructing noncommutative gauge fields starting from some 
symmetry algebra placed at the lattice sites. It might be a quantum symmetry 
algebra, but one can also choose another symmetry algebra. In [39], they choose 
the gauge symmetry algebra. Doing this, it is found that the gauge theory is indeed 
reconstructed if the symmetry algebra is endowed with a co-multiplication. That 
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is to say, that it is extended to a Hopf algebra, essentially. This is a remarkable 
check of the role of the quantum group symmetry in the discretization. 

The previous expression was a discretization of the Stieltjes-Wigert ensemble, 
the discretization of (3) would be: 



(26) Z= J2 e- fl - Y[sinh 2 (u t - Uj ) . 

Ui,...,«AT = — OO i<j 

Note also that, as one can easily see from the observation of the value of the 
coefficients at each point of the lattice, we do not only have a discretization of the 
log-normal or the Gaussian, but there is also a sort of modular transformation -like 

2 2 

in Poisson resummation for example- since the weights are of the type q n = c~ 9s ™ , 
in contrast to the c~" / gs that would come from a direct discretization of the original 
matrix models. This also seems to be analogous to what happens in combinatorial 
quantization, where the symmetry at each site of the lattice is a quantum group 
symmetry. 

The formula (26), coming from the partition function of Chern-Simons on S 3 , 
is easily seen to give the partition function of 2dYM on the cylinder with trivial 
holonomies around the two endpoints of the cylinder [41]. It is also the partition 
function of g-deformed 2dYM on S 2 [8]. As explained in [7, 26], both facts agree 
since these partition functions are the same. Thus, our discretization directly proves 
the equivalence between the Chern-Simons matrix model and 2dYM, including the 
numerical prefactors, which we did not keep track of in (26). This applies to the 
more general cases worked out in [7] and [8] , and in fact it can be generalized to all 
the cases where discretization applies (for example, the Seifert homology spheres 
discussed above). It can also be used in the opposite direction, that is, one starts 
with a discrete model and rewrites it as a continuous matrix model. This can be 
helpful to compute large N limits of such models, since the large N techniques are 
much better understood in the continuous case. We will discuss this in more detail 
elsewhere [38]. 

3. Oscillatory behavior in Chern-Simons matrix models 

3.1. Connection with complex dimensions and fractal behavior. The con- 
nection with complex dimensions [30, 31] can now be easily obtained. We consider 
the Mcllin transform of the g-pcriodic function (15): 

poo 

(27) h(s) = / u(x)x s dx. 

Jo 

Taking into account the following property of Mellin transforms: 

/>oo 

(28) / D(qx)x s dx = q- s h(s), 

Jo 

and considering the g-periodic property (15) : 

27T777 

(29) q- s h(s)=h(s)=>s n = - , neZ. 

log q 

Thus, the Mellin transform of a q periodic function contains infinitely many complex 
poles. From standard results in Mellin asymptotics, this implies that the q periodic 
function may be fractal. This depends on the precise form of the function. Namely, 
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on the residue of the poles. This is a usual situation in Mellin asymptotics with 
complex poles. If the residues decay fast enough with n, then the function will be 
differentiable. 

Another way to arrive at the same conclusion as above, is to notice that 1 + 
?9sin ^27T^|^ in (11) is just one possible deformation among many. The parame- 
ter $ takes its possible values in [—1, 1] just to ensure that the resulting function is 
positive definite everywhere. Of course, the oscillatory term can be easily general- 
ized. For example, one can consider: 

(30) / (x) = c- fe2 >°* 2 *(l + E °» sin ( 27m Iof^) ) ' 

and the integer moments keep their log-normal value. 

3.2. Periodic/Self-similar matrix models. Note that in (30) we are still limit- 
ing the finest scale possible, since the sum stops at n = m and thus the function 
is not a genuine fractal, albeit the oscillatory pattern is certainly much richer than 
in (11). But we can consider a full Fourier series and take m — > oo as long as 
we restrict the coefficients in order to ensure positiveness of the function. This 
restriction is actually very mild. Consider for example, a n — n~ 7 with 7 > 1; then: 

(3!) / W ..^.( 1 + A|±*,(^)). 

is positive definite for a low enough value of the parameter A. One of these functions 
looks as: 




It is plain that for some low enough values of 7, the resulting function is fractal 
[42]. Namely, continuous but differentiable nowhere. An example of such a function 
is plotted in Figure 2. 

Therefore, we have a family of non-polynomial potentials, all of them equivalent 
in the sense that they give rise to the same partition functions. The other possi- 
ble relevant quantities, such as the density of states and correlation functions are 
trivially related [6] . 
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OX) 1 2 3 4 5 



Figure 2. /(x) in (31) with a n = n 1 ,b n = n 2 ,k = 1,A = 
0.5, m = 50 



Going back to the original Chern-Simons matrix model, we are lead to ensembles 
such as, for example: 
(32) 

Z = jj[c-£ + (27rb n Ul ±^ ) ) d "' II ( 2 sinh — 

with suitable a n and b n . At the level of the confining potential we have: 

-Ng s 



( 33 ) V (it) = u 2 - In [ 1 + A ^ a„ sin f l'-/ 

\ n=l ^ 



log<? 



Since V (x) = — Inu) (x) , the derivative of the potential is the logarithmic de- 
rivative of the weight function V (x) = — ^j^j an d thus, non-differentiability of 
lu' (x) implies non-differentiability of V' (x) , since lo (x) is a continuous function. 
Therefore, the potential itself may be a fractal. This possibility depends on the 
choice of the coefficients a n and b n as the previous figures clearly show. A more 
detailed analysis follows from results in Fourier series [42]. While we shall discuss 
this mathematical property elsewhere, recall that if we have a Fourier series such as 
/ ( u ) = Em fl ™ ex P {i- mu ) i where a m have random or pseudorandom phases, then, 
if the power spectrum \a m \ has the asymptotic form: 

(34) \a m \ ~ \m\ as |m| — > oo, where 1 < (3 < 3, 

the graphs of Re/ and Im/ are continuous but non-differentiable, with fractal di- 
mension: 

(35) Df = 1(5-13). 

3.3. Oscillatory behavior of the bare model: density of states. The way os- 
cillatory behavior appears in the matrix models may seem surprising. Nevertheless, 
discrete scale invariance is known to be related to g-deformations with q real [32], 
so its appearance in Chern-Simons matrix models is natural in fact. Furthermore, 
the oscillatory terms seem to be intimately related to an SL(2, Z) invariance. We 
shall discuss this elsewhere, but at least we can readily show now that oscillatory 
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behavior is also a feature even in the original bare model (3). For this, we have just 
to consider the famous expression for the density of states of a Hermitian matrix 
model [1]: 



JV-l 



(36) 



valid for any N. Some features of this quantity have already been studied as they 
are useful in the context of topological strings [4]. However, in this case, one is 
interested in the perturbative expansion in N. This leads one to consider the 't 
Hooft limit where A = Ng s is kept constant while N — > oo. Since this corresponds 
to q ~ 1 , the limit where the g-matrix model approaches a classical one, we lose the 
oscillations in (36), as we shall see. On the other hand, to give a generic expression 
(e.g. the analogous of the semi-circle law for the Hermitian Gaussian model) is 
also not very useful. The problem is that there are apparently no analogues of the 
Planchcrel-Rotach asymptotic formulas [1] and therefore any expression is rather 
cumbersome. But just to get a feeling of the behavior of the density of states 
for arbitrary q, we can carry -employing Mathematica, for example- exact analytic 
computations for various values of N and q, to get an idea of the beautiful behavior 
of the density of states. More precisely, we employ (36) with: 



(37) 



p n { X ) = (-r/v /2+1/4 n(w) 

IE? (i 



1/2 £ 

v=0 



c -k log x with 



, rn=o (i-<? <+1 ) 

That is, the log-normal weight function and the Stieltjes-Wigert polynomials [35]. 
We also go back to the original variables, employing x = e u ~ Ng % and thus, obtain 
several plots of the density of states of the model (3). 




1 .25 
1 

. 75 
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Figure 3. p (u) of model (3), based on (36)-(37) with q = 0.5, TV = 30 

Note how the compact support is much bigger than that of a Gaussian ensemble, 
for example. Actually, it is of the order Ng s , in comparison with the ^/N of a 
Gaussian model. The oscillations that in a Gaussian model are smoothed out, are 
here manifest, as can be seen in Fig. 3 for example. The number of oscillations 
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coincides with N. Crystalline behavior is known to be linked with the presence of 
quantum group symmetries [32] . 




Figure 4. p (u) of model (3), based on (36)-(37) with q = 0.3, N = 10 

With lower values of q the oscillations should be more patent and this is exactly 
what happens in the next plot, Figure 4. Note also how the support grows with 
q. Finally, we go into a higher size where, in contrast to the Wigner-Dyson case, 
the amplitude of the oscillatory behavior remains the same, as can be seen in 
Figure 5. However, note how keeping the same size but taking q towards 1 we are 
approaching the familiar semi-circle like behavior. The oscillations considerably 
decrease in amplitude and so does the support of the distribution. It is therefore, 
as shown in Fig. 6, closer to the Wigner-Dyson case were the oscillations are of 
smaller amplitude and are essentially smoothed out when N is very large. 




. 75 
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Figure 5. p (u) of model (3), based on (36)-(37) with q = 0.5, N = 100 



4. Equivalence of effective superpotentials 



Now we would like to further investigate the possible consequences of the equiv- 
alence that we have found. It seems that this equivalence in Chern-Simons matrix 
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Figure 6. p (u) of model (3), based on (36)-(37) with q = 0.9, N = 100 



models may be related to recent findings in Af = 1 gauge theory [9, 10, 11] (see [43] 
for a review on planar equivalence of related theories). More precisely, these works 
show that different Af — 1 supersymmetric gauge theories can have an effective 
glueball superpotential with exactly the same functional form. In [10] for example, 
it is found that the low-energy effective superpotential of an Af — 1 U(N) gauge 
theory with matter in the adjoint and arbitrary even tree-level superpotential has 
the same functional form as the effective superpotential of a U(N) gauge theory 
with matter in the fundamental and the same tree-level interactions. 

Then, taking into account Dijkgraaf-Vafa [44], the effective superpotential of an 
Af = 1 supersymetric theory is computed by a matrix model whose potential is 
the tree level superpotential of the gauge theory, then the equivalence of effective 
superpotentials, which has been derived employing TV = 1 techniques -like Konishi 
anomaly-, implies the equivalence of the free energies of the corresponding matrix 
models. 

Note that equivalence of free energies is also one of the mathematical results 
obtained in this paper. However, we have discussed Chern-Simons matrix models. 
Do these models have any relationship or implication for Af = 1 theories ? Indeed, 
we have just to recall that at low energies and in four dimensions, IIA theory 
compactified on T*S 3 and with N D6 branes wrapping S 3 reduces to TV — 1 gauge 
theory. Therefore, again due to [44], a Hermitian matrix model describing CS 
theory on S 3 leads to the superpotential of such an Af = 1 theory, which is [45] : 



This role of the Chern-Simons matrix model is already discussed in [18]. Therefore, 
this superpotential can be interpreted as the effective superpotential of an Af = 2 
theory whose tree-level superpotential is the CS matrix model potential, and we 
have just seen that there are infinitely many choices for this potential -of course, 
once again, all of them leading to an identical free energy-, so, all of them lead to 
the same superpotential, namely to (38). 



(38) 




n 
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That is to say, (38) as the effective superpotentiai of an TV = 2 theory whose 
tree-level superpotentiai is: 

(39) V(M) = k 2 \og 2 M + log g(M), 

with a q— periodic g(M) = g(qM), and q = e^ 1 / 2 * 2 = c~ 9s as usual. So, regardless 
the precise form for g (M) , we always end up with a superpotentiai with the same 
functional form. Furthermore, recall that the role of (39) can also be played by a 
discrete potential: 



(40) V(M) = - log J2 c n q^ +n 8(M - cq n ) , 

\n— — oo / 

that, dropping constant terms, comes straightforwardly from (12) and (25). 

Finally, just to mention that in [8], a different Hermitian matrix model for Chern- 
Simons theory was considered: 

oo 2k 

(41) V(M) = M 2 + - a k (~^ S )TrM s TrM 2fe - s . 

fe=l s=0 V s / 
This cumbersome expression, including double-trace terms, comes from the trans- 
formation: 

(42) n ( 2 sinh ^P) 2 = n ^ - f w < 

where / (it) is the logarithm of the second term in (41) [8]. Since this only comes 
from manipulations on the repulsion term of the matrix model, it is independent of 
the weight function. Thus, employing our result, one can consider the generalized 
weight, instead of just the quadratic one and then apply the method in [8], based 
on the trick (42). Therefore, even if one wants to focus on these type of Hermitian 
matrix models, the equivalence result identically holds. That is to say, the addition 
of periodic terms, such as the previously discussed (33), following the M 2 term, 
will not modify the role of (41), as happens with g(M) in (39). 

Thus, the arguments and results of this Section are the same regardless the 
Hermitian matrix models considered. Another different issue is the comparison 
between these Hermitian matrix models, that we shall discuss elsewhere. 

To conclude, note that, while the oscillatory terms may seem at first a bit awk- 
ward, they may be related with the particular form of (38). Namely, with the 
existence of infinitely many domain walls [45]. In addition, these oscillatory terms, 
as mentioned, seem to be related to «SX(2,Z), which is a crucial ingredient in the 
canonical quantization of Chern-Simons theory, and thus in the derivation [5] of 
the original matrix model. 

5. Conclusions and Outlook 

From [6] and the present work, it can be said that numerous Chern-Simons quan- 
tities can be extracted from a certain coarse-grained knowledge of either a Gaussian 
distribution or a log-normal distribution. This coarse-grained knowledge, is noth- 
ing else than a suitable combination of the set of positive integer moments of the 
distribution. This particular combination of the integer moments is the one given 
by the orthogonal polynomials. But it turns out that the matrix models have such a 
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fluctuating weight function that this coarse-grained information -the moments- no 
longer fixes the weight function uniquely, in sharp contrast with polynomial poten- 
tials, for example. Therefore, one has the freedom of modifying the weight function, 
while leaving intact certain amount of information. This reduced information is the 
only one required by the matrix model, and consequently by Chern-Simons theory. 
We have shown the rather large extent to which different superpotentials can be 
equivalent. Mathematically, both the discrete case and the oscillatory -especially 
the ones with very low differentiability- seem interesting on their own. Note that 
the most extreme cases are matrix models with fractal potentials that, nevertheless, 
give rise to physically meaningful objects like Chern-Simons partition functions. 

Regarding possible avenues for further research, discrete matrix models have 
been poorly studied in gauge theory, at least in comparison with continuous models. 
The technology available in these type of models is already quite formidable and, in 
addition, many remarkable combinatorial properties of these models are known [46], 
so their study may prove useful in relationship with concepts relevant to topological 
string theory, like crystal melting. In addition, as discussed above, our methods 
provide a direct understanding of why at the level of the matrix models Chern- 
Simons theory and 2d qYM are equivalent. It also gives a concrete computational 
scheme where one can straightforwardly extend this equivalence between continuous 
and discrete models to other cases, namely, to more general manifolds like Seifert 
homology spheres, and generic torus link invariants. It would be interesting to see 
whether the discrete model that one gets is still reproduced by a BF-type of theory, 
and on which manifold. Hopefully, the equivalence between the continuous and 
the discrete matrix models can be applied to compute large N limits of 2dYM and 
topological strings on more general manifolds as well. Needless to say, our results 
apply also to orthogonal and symplectic groups, so one could consider Calabi-Yau 
orientifolds. 

Interestingly enough, it seems that this non-uniqueness also leads to the equiv- 
alence of effective superpotentials in Af = 1 , through a rather different route from 
[9, 10, 11]. A further understanding of this seems desirable. Another very much 
related question would be to understand if the precise way in which this non- 
uniqueness occurs in the matrix models is related to the specific form of the super- 
potential (38), found in [45]. 
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Appendix A. Notation and Chern-Simons quantities 

We give here some information about the Chern-Simons quantities that appear 
in the text, mainly in (20). For more information, see [5] and references therein. To 
understand the origin of other quantities in (20), one has to take into account the 
constructions of Seifert homology spheres from surgery. Seifert homology spheres 
can be constructed by performing surgery on a link C in S 3 with n + 1 components, 
consisting on n parallel and unlinked unknots together with a single unknot whose 
linking number with each of the other n unknots is one. The surgery data are Pj/qj 
for the unlinked unknots, j = 1, • • • , n, and on the final component, pj is coprime 
to qj for all j = 1, • • • ,n, and the pj's are pairwise coprime. After doing surgery, 
one obtains the Seifert space M = X(^, - • • , This is rational homology sphere 
whose first homology group H\(M, Z) has order \H\, where 

n n 

(43) H = P Y,^> and P =HPi 

3 = 1 Pj 3=1 

Another topological invariant that will enter the computation is the signature of C, 
which turns out to be: 

(44) a(£) = f>ign(£)- S ign0Q. 

For n — 1,2, Seifert homology spheres reduce to lens spaces, and one has that 
L(p,q) = X(q/p). For n = 3, we obtain the Brieskorn homology spheres H{p\,p2,pz) 
(in this case the manifold is independent of qi,q2,qz)- In particular, S(2, 3,5) is 
the Poincarc homology sphere. Finally, the Seifert manifold X(^j, ( m _^)/ 2 ' ^T 2 )' 
with m odd, can be obtained by integer surgery on a (2, m) torus knot with framing 
t. Note that in (20) the weight and root lattices of G are denoted by A w and A r , 
respectively. 

Finally, there is a phase factor in (20) that comes from the framing correction, 
that guarantees that the resulting invariant is in the canonical framing for the 
three-manifold M. Its explicit expression is: 

(H \ n 
- +]Tl2,s( M )--, 
^ / i=i Pl 

where <r(£) is again the signature of the linking matrix of C and s(p, q) is the 
Dedekind sum 

(46) s(p,q) = 1 £oot(^)oot(^). 
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